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£SJ ' We derive a simplified general expression for the two-body scalar Casimir en- 

ergy in generalized separable coordinate systems. We apply this technique to 
the case of radial semi-transparent planes in the annular region between two 
concentric Dirichlet cylinders. This situation is explored both analytically and 
numerically. 



1. Introduction 



In 1948 Casimir 1 predicted that two parallel perfectly reflecting mirrors 
would attract each other with a pressure of P = 7r 2 /240a 4 . Since then 
much work has been done studying a variety of geometries and materials. 
Much of this work has been summarized and referenced in review articles 
by M. Bordag et al 2 and K. A. Milton, 3 and more completely in two books 
by the same authors. 4:5 

This work only concerns itself with the Casimir effect for a massless 
scalar field. In order to proceed we will start with the multiple scattering 
expression for the Casimir energy 



E = ^ / dCTrln(l - G 1 V 1 G 2 V 2 ). (1) 
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Here £ is the imaginary frequency, and is the Green's function referring 
to a single potential Vi . An equivalent expression was first used by Renne 6 
in 1971, and more recently by many others. 7-9 A very good derivation is 
given by Kenneth and Klich. 10 

2. Separation of Variables 

Equation (1) is a fairly complicated formula to work with. We have to per- 
form a 3-dimensional trace of the logarithm of the 1 — G\V\GiVi operator. 
We also have to solve a partial differential equation to find G\ and G 2 . 
However, by working in a coordinate system in which the Hclmholtz equa- 
tion is separable we can greatly simplify this approach. The result will allow 
us to move the trace inside the logarithm, where it will become a simple 
integral, and we will only have to solve an ordinary differential equation to 
find a reduced Green's function for a single coordinate. 

In this section we will find a simplified expression based on a general 
separation of variables using the Stackel determinant. We will follow the 
notation of Morse and Feshbach. 11 

We write the Green's function as a sum of eigenfunctions times a reduced 
Green's function, 



The Mi (£2, £3) is the minor of the Stackel determinant, and the func- 
tions are functions of a single variable related to the scale factors of the 
generalized coordinate system as defined in Morse and Feshbach. 11 The 
X2(^2) and Xzi^s) and a 2 and a 3 are the eigenfunctions and eigenvalues 
determined by the simultaneous set of equations, 



The x eigenfunctions are orthogonal with respect to some weighting func- 
tion p(£2,£ 3 ), 



d£2d£ 3 px2 («2 , a 3 )x2 (a' 2 ,a' 3 )x3 («2 , a 3 )X3 ( (4) 



Using (3), we find that the reduced Green's function in (2) satisfies the 
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differential equation in the single remaining coordinate, 
1 d d 

+ <S> 12 a 2 2 + * 13 o§ + vfafjgitu&C, a 2 , a 3 ) = 5 & (5) 

Working with the Casimir energy written as (1), by expanding the log 
we can write 

oc 

If 00 1 

E = -^ J d(J2-Tr(GiViG 2 V 2 ) s . (6) 

-00 5=1 

The simplification comes if the potentials are functions of only the single 
coordinate £1, with the form Vi(x) — Vi{£,i)/h\. The scale factor hi is 
exactly what is needed to apply the orthogonally condition (4) in performing 
the trace. Finally if the potential consists of two separate non-overlapping 
potentials, we can show 

Tr(G 1 y 1 G 2 y 2 ) s = tr(Si«iS2W2) s - ^ (tr giv ig2 v 2 ) s . (7) 

The interaction Casimir energy can now be written in general separable 
coordinates as 

00 

E= h I d( ^ Ml-tr 5 iiW2). (8) 



3. Casimir Energy for Planes in an Annular Cavity 

As an application we will proceed for the case of two semitransparent radial 
planes in the region between two concentric cylinders, as shown in figure 1. 

This geometry is similar to the wedge geometry first studied in 1978, 12 ' 13 
with a good review by Razmi and Modarresi. However here we include 
circular boundaries in addition to the wedge boundaries. We will enforce 
Dirichlct boundary condition on the inner and outer cylinder. This is similar 
to situations studied by Nesterenko et al 15 ' 16 for global Casimir energies 
for the case of one circular boundary and by Saharian et al 17,ls for the 
local properties of the stress energy tensor for the case of both one and two 
circular boundaries. The radial potentials will be semi-transparent delta- 
function potentials in the angular coordinates, v\(9) = XiS(9) and v 2 {6) = 
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Fig. 1. An annulus with inner radius a, 
outer radius b, and two semitransparcnt 
potentials at 9 = and 9 = a. 



Fig. 2. The contour 7 is defined around 
the positive real line, while not enclosing 
zero. 



\ 2 5{9 — a). This is most similar to the recent work by Brevik et a/, 19 ' 20 and 
Milton et al. 21 

This problem can be solved using separation of variables, leaving £1 as 
the azimuthal coordinate 9. This means we will write our reduced Green's 
function in the azimuthal coordinate, which is different from the traditional 
way of writing the reduced Green's function in terms of the radial coordi- 
nate. From equation (8) we can immediately write 



oo 

£ = ^/dCEMl-tr4 1) ^)- 



(9) 



The Green's function is written in terms of exponential functions that, due 
to the periodicity requirement, give the expression 

m r->\ AiA 2 cosh 2 (77(71- - a)) , 

trg ( - 1 >v 1 g {2) v 2 = — — . (10) 

71 n (2r7sinh?77r + Ai cosh ?77r)(2?7sinh r/7r + A2 cosh ?77r) 

The 77s are the eigenvalues of the modified Bessel equation of purely 
imaginary order, 

d d 



1 1 

-r^r-r— + k r 
or or 



R v {nr) =rj Z R 71 {nr). (11) 



Using the argument principle we can take a complicated sum over eigen- 
values and turn it into a contour integral around the real line as shown in 
figure 2. For this we need a secular function D(rf), which is analytic along 
the real line and has the value zero at the eigenvalues. In this case we define 
R^ko) = then the eigenvalue condition is given by D{rj) — R v (Kb). The 
eigenfunction R v can be written in terms of modified Bessel functions 

R v (nr) = K iri (Ka)Iiri{Kr) - I irj (Ka)K iri (Kr), (12) 
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where we define I v (x) as the part of the modified Bessel function I v (x) even 
in r/. 

The energy per unit length L z can be written as 

00 

d ~ ~ 

— In (K lv (Ka)I lv (Kb) - I lv (na)K irl (Kb)) 



AiA 2 cosh 2 (rjjTT-a)) \ 

(2r) sinh r/n + Ai cosh rjir) (2r] sinh rjir + X 2 cosh rjir) J 

A quick check of this answer is to look at the limit of large inner and outer 
radius, as shown in figure 3. This should then give the answer for a rectan- 
gular piston. For this limit we need the uniform asymptotic expansions of 






K iv and I iv , which are worked out by Dunster. 22 ' 23 We should also redefine 
our dimensionless variables in terms of the dimensionful quantities that will 
appear in the rectangular piston case, rj = rj/a, A = A/a, and d = aa. In 
this asymptotic region we recover the formula for a rectangular piston, 



E 



oc 



^f~^(b - a)) 



— In — 



x In 1 



AiAae- 2 ^ 



• (14) 



(2?y + A 1 )(2^+A 2 ) / 
The contour integral over rj simply ensures that rj 2 = k 2 + (mn/ (b — a)) 2 . 
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4. Numerical Results for Dirichlet Planes 

The Casimir energy in equation (9) is a quickly converging function so it 
should be easy to evaluate. However it can be difficult to evaluate the 77 
eigenvalues, which become functions of the wavenumbcr k and a natural 
number m. We can get around this problem by using (13). We cannot 
integrate along the real line because of the poles introduced when we used 
the argument principle, and we cannot distort the contour to one running 
along the imaginary axis because the integral then becomes divergent. So 
a simple choice is then to let the 77 integration run along the angles of ir/4 

and — 7r/4. Writing tr g^v\g^V2 = A(rj) we have 



e 1 r°° r°° 

Lz 4tt 2 J J 



ReR VJv 9 » ReR VIv +ImR VIv d » ImR V;v + ( TmA(y/iu) 



1^ I 2 \i-R e A(Viv) 

Re R Viv dv Im R sTiv ~ Im R Viv dv Re R sfiv 



2 \ R Vi» 



|2 



xln(l-2ReA(Viv) + \A(Viv)\ 2 ^ |. (15) 



Here we have used the property that R v * = R*, and A(t]*) = A*{rj). The 
value of R^i v (b, k) is obtained as the numerical solution of the differential 
equation. Using this technique we can obtain a numerical energy in about 
1 cpu-second. The results of this calculation are found in figure 4. 

Again we would like to compare to known results, so figure 5 is a graph 
of the ratio of the energies of an annular piston, and a rectangular piston of 
similar dimension. The rectangular piston is constructed so it has the same 
finite width b — a as the annular piston, and the separation distance is the 
mean distance between the annular plates, 

d=— — 2 sin - . (16) 



2 V2, 

The results make a certain amount of physical sense. The energy of the 
annular piston is greater than that of the rectangular piston for small sep- 
aration because the inner edge of the annular piston is closer, and will 
contribute more to the energy. However as the annular piston gets further 
away, the other side of the piston will start to contribute and lower the 
overall energy. In addition we see that the energy for a small piston is much 
closer to that of the rectangular piston for small separations than for a 
larger piston, E ann /E^ cct w 1.004 for b/a = 1.1 vs. E ann /E rcct w 1.23 for 
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Fig. 4. This figure shows the energy per 
length vs the angle between the plates. 
The energy is scaled by the inner radius 
a. 



Fig. 5. This figure shows the ratio of 
the energies of an annular piston to a 
rectangular piston of similar dimension vs 
average separation distance between the 
plates. The separation distance is scaled 
by the finite size of the piston b — a. For 
b/a = 2 only the result for a G [0, it] is 
shown. 



b/a = 2. In both cases the value approached in the plateau in figure 5 is 
very close to the ratio of the energies of a flat plate to that of a tilted plate 
predicted by using the proximity force approximation. 
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